Response to Comment on 'Undamped electrostatic plasma waves' [Phys. Plasmas 19, 

092103 (2012)] 
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Numerical and experimental evidence is given for the occurrence of the plateau states and con- 
comitant corner modes proposed in Q. It is argued that these states provide a better description 
of reality for small amplitude off-dispersion disturbances than the conventional Bernstein- Greene- 
Kruskal or cnoidal states such as those proposed in 0]. 

PACS numbers: 52.20.-j; 52.25.Dg; 52.65.-y; 52.65.Ff 



Since the publication of the original Bernstein-Greene- 
Kruskal (BGK) paper [l|, which described ways to con- 
struct a large class of nonlinear wave states, there has 
been an enormous literature that speculates about which 
of these states might occur in experiments, in nature, 
and in numerical simulations in various situations 0] ■ In 
his Comment [3| , Schamel has written a broad spectrum 
diatribe touching on many points; we agree with some of 
these points, in fact, some of them were originally advo- 
cated by some of us. In particular, we are in agreement 
that off-dispersion excitations, made possible by particle 
trapping, are important. However, we will confine our 
response to those of his comments that are relevant to 
our paper, which is about the construction of an appro- 
priate linear theory that describes the small amplitude 
limit. Schamcl's claim is that our lack of trapped par- 
ticles invalidates our analysis even in this limit. Rather 
than reiterate the points made in our paper, we rebut his 
claim by providing further numerical and experimental 
evidence. 

Our point is best made by a glance at Fig. @] below 
that shows late time simulation results for two runs: 
an off-dispersion case (Run B), with the wide nearly x- 
independent plateau that we proposed for the distribu- 
tion function, and an on-dispersion case (Run A) , depict- 
ing the more conventional BGK or cnoidal wave state. 
Evidently, the more conventional BGK or cnoidal wave 
state is not the best description of the late time dynamics 
of the off-dispersion simulation. 

In the following we will describe our nonlinear simu- 
lations that produced Fig. 2] and give comparison to our 
theory of j 1) . This will be followed by a discussion of 
some experimental results that further provide evidence 
for the validity of our theory and, in particular, the ex- 
istence of corner modes. 

As in Ref. [H, we make use of a Eulerian code 0-0] 
that solves the Vlasov-Poisson equations for one spatial 
and one velocity dimension: 

2L + v ^l- E ^l = o ?E = i-[fdv (i) 

dt dx dv dx J 



0N-DISPERI0N MODE 




4000 



4000 



FIG. 1: (Color online) Time evolution of the fundamen- 
tal spectral component of the electric field Ek(t) of an on- 
dispersion mode (Run A at the top) and of an off-dispersion 
mode (Run B at the bottom). The vertical dashed lines rep- 
resent t a g. 



where / = f(x,v,t) is the electron distribution function 
and E — E(x,t) the electric field. In ([1]), the ions are a 
neutralizing background of constant density hq = 1, time 
is scaled by the inverse electron plasma frequency ojp 1 , 
velocities by the electron thermal speed v t h , and lengths 
by the electron Debye length \rj. For simplicity, all the 
physical quantities will be expressed in these characteris- 
tic units. The phase space domain for the simulations is 
T> = [0, L] x [— u a , i> max ]- Periodic boundary conditions 
in x are assumed, while the electron velocity distribution 
is set equal to zero for |u| > w max = 6. The x-direction 
is discretized with N x = 256 grid points, while the v- 
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FIG. 2: (Color online) Resonance peak for Run A (top) and 
Run B (bottom); the red- vertical lines indicate the value of 
the driver phase velocity vj, for Run A and B, respectively. 



direction with N v = 12001. 

In our previous simulations of Ref . [H , the initial equi- 
librium consisted of a velocity distribution function with 
a small plateau; however, here we assume a plasma with 
an initial Maxwellian velocity distribution and homoge- 
neous density. We then use an external driver electric 
field that can dynamically trap resonant electrons and 
create a plateau in the velocity distribution. This is 
the same approach used in the numerical simulations of 
Ref. [1-G3 and in the experiments with nonneutral plas- 
mas in Ref. fli"l |. 

The explicit form of the external field is 

E D (x, t) = g{t)E DM sin (kx - u D t) , (2) 

where E DM is the maximum driver amplitude, k = 2tt/L 
is the drive wavenumber with L the maximum wave- 
length that fits in the simulation box, uj d = kv^^ is the 
drive frequency with v<p D the driver phase velocity, and 
g(t) = [1 + (t — t)"/At"] _1 is a profile that determines 
the ramping up and ramping down of the drive. The 
external electric field is applied directly to the electrons 
by adding E D to E in the Vlasov equation. An abrupt 
turn-on or turn-off of the drive field would excite Lang- 
muir (LAN) waves and complicate the results. Thus, we 
choose n = 10 so g(t) amounts to a nearly adiabatic turn- 
on and turn-off. The driver amplitude remains near E DM 
for a time interval of order At centered at t = r and is 
zero for t > t Q g ~ r + At/2. We will analyze the plasma 
response for many wave periods after the driver has been 
turned off. 

We simulate both the excitation of an on-dispersion 



FIG. 3: (Color online) Phase space contour plot of the elec- 
tron distribution function at t = t max for Run A (top) and 
Run B (bottom). Black dashed and black dot-dashed lines in- 
dicate the driver phase speed v<f, D and the mode phase speed 
v$, respectively. 



mode (Run A), a mode for which (k, v^) is on the thumb 
curve of Fig. 1 of Ref. and an off-dispersion mode 
(Run B), that is off the thumb curve. The excitation of 
the on-dispersion mode is obtained through an external 
driver with k = 7r/10 and v ( j >D — 1.45, while for the off- 
dispersion mode we set k = 0.7 and = 1.5. The 
maximum driver amplitude E DM =0.01 has been chosen 
for each simulation in such a way that At ~ 10r t , with 
Tt being the trapping period fl2l ]. Finally, the maximum 
time of the simulation is t max = 4000, while the driver is 
zero for t > t oS = 2000. 

Figure [1] shows the time evolution of the fundamen- 
tal spectral component of the electric field, Ek(t), for 
Run A (top) and Run B (bottom). In both plots one 
can see that after the driver has been turned off at t a ff 
(indicated by the red-dashed lines in the figures), the 
electric field oscillates at a nearly constant amplitude. 
Figure depicts the resonance peaks for Run A (top) 
and Run B (bottom), obtained through Fourier analysis 
of the numerical electric signals performed in the time in- 
terval t f[ < t < t max , i.e., in the absence of the external 
driver. From these two plots it is readily seen that the on- 
dispersion mode propagates with phase speed very close 
to the driver phase velocity v^ d , whose value is indicated 
by the vertical red-dashed lines, while the phase speed of 
the off-dispersion mode is shifted towards a lower value 
with respect to v^^ . This shift is predicted by our theory 
of Ref. [H and these results provide qualitative evidence 
for its validity; subsequently, we will show quantitative 
agreement. 



3 




ON-DISPERSION MODE 




FIG. 4: (Color online) Surface plot of the electron distribution 
function at t = t max for Run A (top) and Run B (bottom) . 



The main differences between the on-dispersion and 
off-dispersion modes can be appreciated by examination 
of the phase space contour plots of the electron distri- 
bution function shown at t — t max in Fig. |3] For Run 
A (top) a well defined trapping region that propagates 
in the positive x-direction is visible. The black-dashed 
and black dot-dashed lines in the figure represent the 
phase speed of the driver v<j >D and excited mode v^, re- 
spectively. For Run A it is easily seen that and 
are almost identical, meaning that the region of trapped 
particles generated by the external driver survives even 
when the driver is off and that this region streams with 
a mean velocity close to that of the driver. The physi- 
cal scenario appears quite different for the off-dispersion 
mode of Run B (bottom) . Here we observe a rather wide 
nearly x-independent plateau (the orange region of the 
plot) that is substantially wider than the separatrix for 
the trapped particles (small dark region at velocity close 
to v ~ 1.4). Moreover, here the values of v$ D and v^, are 
well-separated, meaning that the excited mode oscillates 
with a frequency smaller than that of the external driver. 

The differences between Run A and Run B can be fur- 
ther appreciated by looking at the electron distribution 
function surface plots of Fig. @] For Run A (top) we 
observe a trapped region modulated in the spatial di- 
rection, while for Run B (bottom) we see a flat region 
whose velocity width appears to be independent of x. 
The on-dispersion Run A resembles the more conven- 
tional BGK type solution like Schamel's, while the nearly 
x-independent off-dispersion plateau of Run B is very dif- 
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FIG. 5: Velocity dependence of f(L/2,v) at t — t max for 
Run A (top) and Run B (bottom); black dashed and black 
dot-dashed lines indicate and v^, respectively 



ferent, it being more like a quasilinear plateau. Since the 
on-dispersion case has no frequency shift, it appears that 
the trapping dynamics is dominated by a single wave 
and a BGK type solution is to be expected. For the off- 
dispersion case, where there is a frequency shift between 
the driver and the ringing wave, the trapping dynamics 
may involve multiple waves with different phase veloci- 
ties. The interaction between these waves could be caus- 
ing a band of chaotic dynamics that re-arranges phase 
space to provide a more quasilinear type of plateau. 

Quantitative evidence for our theory can be extracted 
from Fig. [5j which shows /(x — L/2, v) as a function of 
v for Run A (top) and Run B (bottom). Again, black 
dashed and black dot-dashed lines indicate and v^, 
respectively. Also here the phase velocity shift for the 
off-dispersion mode is evident; by taking into account 
the uncertainty due the finite time resolution of the sim- 
ulations, we can estimate the interval in which the value 
Av^ um ^ := — v^ d of the phase speed shift falls. This 

gives -0.095 < Av^ um) < -0.0104. Using this we can 
compare the phase velocity shift obtained for Run B to 
the analytical prediction using the "rule of thumb" of 
Eq. (23) in Ref. the theoretical expectation for the 
phase velocity shift of the off-dispersion mode of Run B 
is Av { * h) ~ -0.0946 (with a value Av$ h) ~ -0.0933 ob- 
tained by increasing the resolution by two order of magni- 
tude in velocity), in very good agreement with the value 
obtained from the simulation. Thus our theory not only 
predicts the qualitative direction of the phase velocity 
shift, it gives a very good quantitative value. 

To conclude, we show that there already exists pub- 
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lished qualitative experimental evidence [ll| for the valid- 
ity of the theory we proposed in Q . Figure 6 of Ref. [ll| 
depicts the plasma response to a drive at a spread of 
frequencies. In Fig. 6(b) the larger peak on the right cor- 
responds to the Trivelpiece- Gould mode, which in this 
experiment corresponds to the LAN mode of the thumb 
curve, while the small peak on the right corresponds to 
the Electron-Acoustic wave. Thus, frequencies between 
these peaks and above the LAN peak correspond to off- 
dispersion modes. In addition, at the bottom of Fig. 6(b) 
is given an indication of the frequency shift between the 
plasma response, corresponding to a ringing mode, and 
that of the drive. Since k is fixed, this frequency shift is 
equivalent to a shift in the phase velocity. Observe that 
the frequency shift is positive within the thumb curve 
(between the peaks) and is negative above the LAN peak 
The directions of these shifts can be inferred from our 
theory, cf. the rule of thumb, Eq. (23) of Ref. [![. In this 
equation k 2 = M gives a frequency on the thumb curve, 
while a frequency of a corner mode, an off-dispersion ex- 
citation, is obtained by seeking a root with the addi- 
tion of the plateau contribution (the remaining term of 
Eq. (23)). Because the rule of thumb gives the local shape 
of the plateau contribution, it is not difficult to infer the 
direction of the frequency shift relative to the drive. A 
straight bit of reasoning using the rule of thumb implies 
frequencies within the thumb curve and those above the 
LAN mode should shift in precisely the directions seen 



in the experiments. 

More details about the simulations discussed here and 
further experimental verification of our theory will be the 
subjects of future works that are presently under prepa- 
ration. 

Postscript: In his response to the first version of the 
present Response, Schamel modified his Comment |3/ in 
an attempt to use our numerical and experimental re- 
sults to substantiate his case. We are not convinced by 
his arguments and stand by our original conclusion of 
Q/; viz., corner modes, under the circumstances we de- 
scribed, provide a better description of computational and 
experimental results than cnoidal/BGK modes. 
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